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Q ■ Abstract 

CN I A proof of the Jordan canonical form, suitable for a first course in linear 

algebra, is given. The proof includes the uniqueness of the number and sizes of 
the Jordan blocks. The value of the customary procedure for finding the block 
I generators is also questioned. 
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The Jordan form of linear transformations is an exceeding useful result in all theo- 
^ . retical considerations regarding conjugacy classes of matrices, nilpotent orbits and the 

Jacobson-Morozov theorem. A classical reference for this topic is Smirnov's book [7, 



m 
(N 



p. 245- 254]. There is a very well known proof due to Fillipov [3], which is also given in 
Strang's book [8, p. 422-425]. The American Mathematical Monthly has published at 
least six proofs of the Jordan form over the years: [1, 2, 4, 5, 6, 9]. The justification 
for approaching the subject yet another time can only be the clarity and brevity of 
the presentation and a new criterion for the uniqueness of the number and sizes of the 
Jordan blocks. This note gives such a proof, which has the added advantage that the 
most important parts can be taught in a first course on linear algebra, as soon as basic 
ideas have been introduced and the invariance of dimensions has been established. It 
is thus also a contribution to the teaching of these ideas. 

Although extensive work has been done in [10] regarding this circle of ideas, the 
method given in this note provides a very simple algorithm whose efficiency is shown 
through worked examples. 
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In view of the algorithm given in this note, and the examples given below, it is not 
clear to us why a precise determination of the block generators is needed, although, for 
the sake of completeness, we have discussed this aspect too- at the expense of increasing 
the level of exposition. 

It would be very desirable to compare the computational complexity of computing 
the Jordan canonical form, using the algorithm given in this paper, with other algo- 
rithms, which programmes like Maple and Mathematica use to determine the Jordan 
form. 

As is well known, the main technical step in establishing the Jordan canonical form 
is to prove its existence and uniqueness for nilpotent transformations. We will return 
to the general case towards the end of this note. 

Let A he a nilpotent transformation on a finite dimensional vector space V, let v 
be a nonzero vector in V and n the smallest integer such that A^v = 0. 
Proposition 1 The vectors {A^v : < i < n} are linearly independent. 
Proof. Take an expression 

n-l 

J^CiA'v^O, (*) 

1=0 

in which the number of non-zero coefficients is as small as possible. If the coefficient 
Cj is the non-zero coefficient of largest index j, then multiplying by A""^, we obtain 
an expression like (*) of smaller length. So in (*) every q with i < j is 0. Therefore 
CjA^v — and therefore A^v — 0, with j < n — 1, which contradicts the choice of n. 
This proves the claim, ■ 

Proposition 2 Let R{A) be the range space of A and N{A) be the null space of A. Let 
{A{vi) : i — 1, . . . ,r} be a basis of the range space. Let {uj : j = 1, . . . , s} be a basis of 
the null space of A. Then {vi : i — 1, . . . ,r,nj : j — 1, . . . , s} is a basis of the vector 
space V. 

r r 

Proof. Let v & V. So A{v) = ^ CiA{vi). Therefore v—^ CiVi belongs to the null space 

1=1 i=l 

of A, hence it is a hnear combination of the {vi} and {uj}. To see that these vectors 
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are linearly independent, suppose ^ j + ^ djUj = 0. This gives ^ cA{vi) = and 

i=l j=l i=l 

by linear independence of the vectors A{vi), we get Ci — 0, i — 1, . . . ,r. The linear 
independence of {rij} then shows that dj = 0, j = 1, . . . , s. m 
Proposition 3 V is a direct sum of cyclic subspaces. 

Proof. We prove this, as in the standard proofs [7,8], by induction on dimension. The 
null space of ^4 is a non-zero subspace and therefore the range space of ^4 is a proper 
subspace of V. If this is the zero subspace, then a basis of V gives the decomposition 
into cyclic subspaces. So suppose that R{A) is a nonzero subspace. It is an A invariant 
subspace. By induction on dimensions, it is a direct sum of cyclic subspaces, with 
generators Vi, i = 1, . . . , k, and basis A^Vi, < j < rii, and j = 0. Let Vi = Awi. 

So A^Vi — AA^Wi shows, using Proposition 1, that the vectors A^Wi, < j < rii are 
linearly independent. Also A^^'^^Vi — A^^'^'^Wi — 0, so A^^'^^Wi — A^^Vi belong to the 
null space of A. 

By Proposition 2, if we enlarge A^^Vi, i — 1, . . . ,k, to b. basis of the null space of A 
by adjoining independent vectors rii, . . . ,ni in the null space of A, then A^Wi,0 < j < 
ni,0 <i<k, A"'*Vi,i — 1, . . . ,k, rii, . . . ,ni form a basis of V. 

Therefore, the cyclic subspaces generated by Wi, i = 1, . . . , k and the one- dimensional 
subspaces generated by n^, 1 < r < I give a direct sum decomposition of V into cyclic 
subspaces. ■ 

Prom this description, it is clear that in each summand only A"'^'^^Wi — A^^Vi 
contributes to the null space of A in that summand and therefore the number of 
summands in the above given decomposition is the dimension of the null space of A. 
Corollary Let di — dim {N{A\R{A^)), i = 0, 1, . . . , n), where n is the smallest positive 
integer so that A^ = 0. The differences do — di,di — d2, ■ ■ ■ , c?n-i — dn give the number 
of Jordan blocks of sizes 1, 2, . . . , n. 

Proof. As shown in the proof of Proposition 3, the number of summands in the Jordan 
decomposition is the dimension of the null space of A. Therefore the number of blocks 
of size> 1 is dim(A^(A)). Applying A removes all blocks, if any, of size 1, and so the 



number of blocks of size> 2 is dim{N{A\R{A))) = di. Continuing, we get that di is 
the number of blocks of size> i, i — 1, . . . ,n. Therefore the difference cij_i — d^ gives 
the number of blocks of size i, for i — 1, . . . ,n. m 
Examples 



1. Let A be any nilpotent upper triangular matrix whose entries to the right of 
the main diagonal give a non-singular matrix. Then the null space of ^4 is 1 
dimensional and therefore the canonical form of A consists of only one block. 

In particular, the matrices 

2 
1 
-1 

-2 


12 3 4 

7 6 5 
8 9 
10 


are conjugate matrices as they are conjugate to 

1 
1 
1 
1 




and 



2. If 



10 1 

10 

1 





then N(A)works out to be 1 dimensional, so there is only 1 Jordan block. 



Also 



A' 



1 









so 



N{A^) 



X 

y 

z 




and, as = 0, a basis of N{A'^)/N{A^) is 

" " 



• 
1 

Therefore, this must be a generator of the block. 



3. Let 



A 



2 2 1 

2 11 

2 2 

4 



The eigenvalue 2 is of multiplicity 3, so the generalized eigenspace V(2) is 3— 
dimensional, whose basis works out to be (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0) and the 
matrix of vl|V(2) is therefore 

"202 
2 1 
2 



We have 



(^-2/)|% 



2 
1 




Let A = {A-2I)\V^2)- This gives do = dim(Ar(i) = 2, di = dim{N{A\R{A)) = 1, 
d2 = dim{N{A\R{A^)) = 0. 

Therefore A has do — di = 1 block of size 1 and di — d2 = I block of size 2. 
The Jordan form of A is therefore 




1 




and of ^|V(2) is 

"200 
2 1 
2 

The eigenspace for eigenvalue 4 is one-dimensional. Therefore, the Jordan form 
of A is 

" 2 
2 10 
2 
4 
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In view of such examples, it is not clear to us why an algorithmic procedure is 
needed to find the precise generators of the various blocks, because all one needs to 
find the form of the Jordan blocks is to compute the invariants dj. Nevertheless, for 
the sake of completeness, we outline such a procedure- at the expense of increase in 
level of exposition. 
Step 1: 

Find all eigenvalues. For an eigenvalue A, compute the generalized eigenspace 
corresponding to A. Although, one needs to compute only all vectors annihilated 
by {A — A/)^'™^, it is algorithmically better to compute the vectors annihilated by 
{A — A/)" = 0, where n is the multiplicity of the eigenvalue A in the characteristic 
polynomial of ^4. So, by working in the generalized eigenspace for A, and replacing 
(^4 — XI) by A, we may assume that ^4 is a nilpotent transformation of index < n. 

From now on, we assume that A is a. nilpotent transformation defined on a vector 
space V 
Step 2: 

Find the number and sizes of blocks of this nilpotent transformation according to 
the algorithm given below: it is a restatement of the Corollary given on p. 3. This is 
the most important step, which is needed to complete the next step efiiciently. 
Algorithm for finding the Jordan Form 

For a nilpotent transformation ^4 on a finite dimensional vector space V, let N be the 
smallest integer such that A^ = 0. Let di = dim {N{A\R{A')), i = 0, 1, . . . , A^). 
The differences 

^^0 ~ di, di — d2, ■ ■ ■ , d]sf-i — d]sf 

give the number of Jordan blocks of sizes 1, 2, . . . , A^. 

Step 3: Algorithm for finding the block generators 

Call a nonzero vector v is of height n if n is the smallest integer so that A^^v) = 0. 

The vector space spanned by v, Av, . . . , A^'^v is n-dimensional. A block of size n is 

an A-invariant subspace generated by a vector of height n. 
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Let n be the size of the largest block. Choose a basis of N {A^) / N {A^~^) . This is a 
non-zero space, because there exist blocks of size n. The smallest 74-invariant subspace 
of the preimages gives a direct sum of blocks, each of size n. Call this space W^. 

Let m be the size of the block immediately below n. Consider N [A"^) / N [A"^~^) . 
Find a basis of N{A'^\W{)/N{A"'-^\W{). 

Let wi, . . . , be the preimages of these basis elements; they are all of height m. 

Extend this basis of N {A'^\Wi) / N {A'^-^\Wi) to a basis of N {A"^) / N {A'^-'^) by 
adjoining independent elements with preimages Vi, . . . , v^. 

The smallest >l-invariant subspace spanned hj vi, . . . ,Vs - call it W2 has intersec- 
tion with W\. 

Let H/^i © FF2 = Let I be the size of the blocks, if any, just below m. Extend 
a basis of N{A^\W^)/N{A^-^\W^) to a basis of N {A^) / N {A^'^) . As before, we will get 
the required number of blocks of size / complementary to Wi © W2 = W3. Continuing, 
this will give a Jordan decomposition. 
Explanation 

Step 3 is based on the following observations 

1. If 1^ is a direct sum of blocks and the size of the smallest block is n and < j < n, 
then the null-space of A^ in W is the range space of A"~^ in W. 

2. If is a direct sum of blocks of size n, generated by vectors Vi, . . . ,Vk - all of 
height n, then these vectors are linearly independent in N {A"') / N {A^~^) . 

Conversely, if vectors wi,...,wi are in N{A^) and their images in the quotient N{A^) /N{A'^~^) 
are linearly independent, then the smallest 74-invariant subspace generated by Wi, ... ,wi 
is a direct sum of blocks of size n, with generators wi, . . . ,wi. 
Example: 
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Using the above algorithm, the reader can check that if 
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then the smallest integer n so that — is 6. There are two blocks, of sizes 1 and 6 
respectively, generated by 
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and 
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A final remark on applications: A main application of the Jordan form in differen- 
tial equations is in computation of matrix exponentials. However, it is computationally 
more efficient to calculate the matrix of A relative to a basis of generahzed eigenvectors- 
not necessarily given by cyclic vectors -and compute its exponential relative to this ba- 
sis; finally, conjugating by the change of basis matrix gives the exponential of A. 
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